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Abstract. We introduce a general scheme to realize perfect storage of quantum 
information in systcms of interacting qubits. This novel approach is based on 
global external controls of the Hamiltonian, that yield time-periodic inversions in 
the dynamical evolution, allowing a perfect periòdic quantum state recontruction. 
We illustrate the method in the particularly interesting and simple case of spin 
systems affected by XY residual interactions with or without static imperfcctions. 
The global control is achieved by step time-inversions of an ovcrall topological 
phase of the Aharonov-Bohm type. Such a scheme holds both at finite size and 
in the thermodynamic limit, thus enabling the massive storage of arbitrarily large 
numbers of local states, and is stable against several realistic sources of noise and 
imperfections. 
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1. Introduction 

In the attempt to pave the way to the realization of scalable schemes for quantum 
computation, much theoretical work has been recently aimed at developing suitable 
strategies for the efhcicnt processing and the coherent transfer of quantum information 
0211210] 03- Besides these two fundamental aspects, a further crucial requirement 
for the realization of scalable quantum computers is the possibility to store quantum 
data on time scales at least comparable to those needed for the computational process. 
In particular, it is very important to introduce systems acting as stable and robust 
quantum memòries that recover and conserve large sets of quantum states that would 
be otherwise usually lost in very short times, due to quantum diffusion and decoherence 



To ensure stable information storage in a quantum register, many different noise- 
evading schemes have been proposed |3 El d El El EÏÏl HH [T21 ESI H3| - All these 
works can be roughly classified in two different groups: the first includes schemes 
based on some error correction technique. The remaining proposals exploit somc 
intrinsic property of the quantum register that leaves some specific subsets of quantum 
states unaffected along the temporal evolution (decoherence- free subspace schemes). 
These latter approaches provide, in principle, the complete solution to the problem 
of quantum information storage, but, unfortunatcly, they are extremely sensitive to 
almost any source of imperfection. On the contrary, the schemes based on quantum 
error correction techniques are characterized by a dynamics that allows at any time 
the unambiguous reconstruction of the initial information. The main trouble with 
error correction techniques lies in the fact that usually only very few states can be 
effectively accessed to store information and, hence, relatively large arrays of qubits 
are needed to memorize relatively small amounts of information. 

In the present paper we introduce a new approach to quantum state storage 
based on the idea of time-controlled periòdic dynamics that allows a perfect, periòdic 
reconstruction of a genèric initial state. To demonstrate and describe it, let us first 
recali that an ideal quantum register can be considered as a set of isolated identical 
qubits subject to a local Hamiltonian 



where B is the half gap between the two energy levels of each spin^S]. However, 
in realistic situations, the register is subject to noise caused by disorder in the local 
gap and by interactions both with the substrate environment and between the qubits. 
Typically, at least at sufRciently low temperatures or sufficiently weak coupling with 
the background substrate, the interaction with the external environment takes place 
on time scales much slower than those associated to the computational process. For 
instance, in the case of quantum gate operations with hypcrfmc levels in trapped 
ions, the ratio t gat e/tdecoh can be as small as 10~ 9 03 El- Then, the corrupting 
effects that take place on the same time scales of the computational processes, 
and thus need to be addressed first, are those due to the unavoidable presence of 
the residual, deterministic and/or random inter-qubit interactions 03 E|. These 
interactions cause, in general, fast quantum diffusion and, as a consequence, the 
complete corruption of the information one wishes to store and process. Therefore, 
even before considering the effects of thermal and environmental decoherence, realistic 
quantum registers in the presence of irreducible noise and imperfections must be 








(1) 
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described by Hamiltonians of the form |151 IT%| : 

Htot = Ho + H err , (2) 

where H err is the residual inter-qubit Hamiltonian. In several physical situations, it 
can be described by XY interaction tcrms: 

Herr = ~ + Vi) Ut a ï+l + H - C ~) ' ( 3 ) 

i 

where the random variables of vanishing mean r\i are the local imperfections in the 
global, averaged nearest-neighbor coupling amplitude A, and er^ = a% ± icr v k . Such a 
site-dependent XY model applies immediately to spin-1/2 based quantum registers 
(such as in NMR devices) , but gives as well an effective description of a register based 
on hopping and/or interacting partides on a lattice, in the presence of an energy gap 
such that only two local states on each lattice site can be considered Jü] , so that the 
material partides are mapped in spin-1/2 systems. The storing scheme that we shall 
describe in the present work holds in general for many classes of inter-qubit interaction 
Hamiltonians, for instance XXZ interactions, but in a more limited range of validity 
when dealing with imperfections 20 . Therefore, in the following, we will restrict our 
analysis to physical systems and regimes such that the residual inter-qubit couplings 
can be well described by XY interaction Hamiltonians, including local imperfections. 

As already mentioned, the presence even if A <C B and rjt = V i, 

rapidly destroys the storing capacity of the register: immediately after having stored 
an initial state \ip(t — 0)) = |%), such state starts to evolve and diffuse indcfmitcly. To 
illustrate and measure this effect, we can follow the behaviour of the time-dependent 
two-dimensional overlap Td{t) = KxIV'W)! 2 f° r different stored initial states and 
physical situations |21| . Let us then consider a closed (periòdic boundary conditions), 
unmodulated chain with XY nearest neighbour residual interactions between the 
spins, and let us introduce the set of one-magnon, locally excited states 

i**>=ií«i>nu'>' ( 4 ) 

where the product involves all the qubits of the register except the one at site d. In 
Fig. 1 we take the initial state \x) to be l^o^ i- e - the one-magnon excitation placed 
on site d = 0, while in Fig. 2 the initial state is taken to be the linear combination 
with equal weights of one-magnon excitations placed on sites d = 1 and d = — 1. 
By looking both at Fig. Q and Fig. J2J), we see that, irrespective of the different set 
of physical parameter and stored states one considers, the quantum state diffusion 
grows indefinitely in time and the initial state is never recovered at any site of the 
lattice. Consequently, the one-dimensional hdelity of the initially stored state at a 
fixed position d in the lattice quickly decrease, undergoing partial revivals, as shown 
in Fig. ©. It is worth noticing that in a non modulated register the hdelity of an 
initial superposition is more robust against quantum diffusion. This intrinsic property 
is helpful in the case of modulated (controlled) registers. 

2. Storing quantum information by topologically modulated spin chains 
and time-inverted quantum dynamics 

To overcome this major obstacle in the construction of prima facie working models of 
a quantum register we will introduce a global (that is, site-independent) control of the 
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Figure 1. Unmodulatcd rcgister: contour plot showing the evolution of the 
bidimensional overlap T d (t), starting from the initial state l^o)! as a function 
of the distance d from site (j/-axis) and of the dimensionless time Aí in atòmic 
units fi = 1 (x-axis). Hcrc B = 100A. The value of the overlap increase from 
white {F d (i) = 0) to black (JF d (t) = 1). 




Figure 2. Unmodulated register: contour plot showing the evolution of the two- 
dimensional overlap for the initial state |x) = (1*1) + |^ — i))/v2, as a 

function of the distance d from site (y-axis) and of the dimensionless time Aí 
(x-axis) in atòmic unit h = 1. The value of the overlap increases from white 
ÇF d {t) = 0) to black (JF d (t) = 0.5). 



total Hamiltonian H tot able to force a periòdic time reconstruction of the initial state. 
In the presence of a genèric time-dependent Hamiltonian, evolution of a quantum 
state can be determined in the most general terms by resorting to the Dyson series 
representation. Obviously, the structural complexity of the Dyson evolution integral 
does not allow to identify all the possible Hamiltonian dynamical evolutions yielding 
perfect, time-periodic quantum state reconstruction. However, the Dyson series can 
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Figure 3. Unmodulated register: One-dimensional fidelity To{t) at fixed latticc 
site d = 0, as a function of the dimensionless time At for difïerent initially stored 
states. The sòlid line represents the evolution of the fidelity with initial state taken 
to be the single-magnon excitation |^o} localized at site d = 0. The dashed line 
represents the evolution of the one-dimensional fidelity for the initial superposition 
state |x) = (|\£'i) + | , I' — 1))/%/2. Notice that superpositions are more robust against 
quantum diffusion than non-superpositional, localized excitations. 



be easily resummed when the Hamiltonian enjoys the property to self-commute at 
difïerent times: [H(t), H(t')] = for all t and t' ^ t. Thus, our first prescription 
concerns a quantum register described by a time-periodic modulated Hamiltonian, 
commuting at all times, either when applied to all possible states, or for a subspace 
of the whole Hilbert space of states. 

Under the time-commutation hypothesis, we can determine a complete set of 
states {\a}} that are eigenstates of the time-dependent Hamiltonian at all times even 
if the corresponding energy eigenvalues are time-dependent functions £ a (t). Because 
the energy eigenstates form a complete basis set, the initial state \x) can be written 
as a linear combination: \x) — ^2 a c a \a), with c Q = (a|x), and the sum runs over the 
complete set of eigenstates. It is then simple to write the evolution of the initial state 
after a time t > 0, 

m)) = Ç c a cxp (-i jí e a (r)d^j \a) . (5) 

From Eq. (J5J it is immediate to see that if at a certain time T > all the integrals 
in the sum are equal or differ from each other by integer múltiples of 27r, then the 
initial state \x) is perfectly reconstructed, but for an irrelevant global phase factor. 
This is then the second requirement that one needs to impose on the time-modulated 
dynamics in order to realize perfect time-periodic quantum state storage. 

These two bàsic requirements can be implemented successfully for spin systems 
with XY interactions by a scheme of quantum control based on time-inverted 
dynamical evolutions realized by suitably engineered time-dependent global phase 
factors. We then consider the following situation. In the presence of a simple geometry 
realized by placing a tiny solenoid in the center of a circular ring (periòdic boundary 
conditions) of spins sitting at regularly spaced lattice sites, all the nearest-neighbor 
interaction amplitudes become complex by acquiring the same site-independent, global 
phase 9 proportional to the magnètic flux linked to the ring: 9 cx (p/N, where N is 
the total number of qubits in the ring. The solenoid field can then be modulated in 
time to achieve the control needed for perfect quantum state storage. In the presence 
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of a time-variable linked flux, the total Hamiltonian Eq. Q becomes time-dependent 
and reads: 

H tot (t) = -5> + 7?i ) (e <fl <V*ín + H - C ) + fl E CT i ■ ( 6 ) 

z i 

Incidentally, we note that this scheme can apply as well to more general systems 
of hopping material partides, thanks to the "Peierls" effect, i.e. the fact that if 
charged partides are in the presence of a linked magnètic flux, the real-valued hopping 
amplitude between partides (in spin language, the nearest-neighbor coupling), is 
transformed in a complex- valued quantity |22| . 

Clearly, not all time modulations of the phase factor can realize the dcsircd pcrfcct 
state storage. The first constraint to be imposed is commutativity at different times: 
[H to t{t),H tot {t')] = 0. This property is verified if and only if 9{t) - 9(t') = kir, 
with k integer. This implies that during the entire time evolution the phase must 
be modulated in regular periòdic jumps (steps) between two constant vàlues 9q and 
6{) + 7r (step-phase modulation). 

Concerning the property that at a certain given time T, all the time integrals 
appearing in the expression Eq. © must be equal (or differ by a trivial phase factor 
integer múltiple of 2ir), let us observe that, independently of the vàlues of 9, the local 
term Hq (ideal register) commutes with the XY residual interaction terms in the time- 
dependent Hamiltonian Eq. JBJ. Hence, there exists a complete set of eigenstates of 
H c (t) that are as well simultaneous eigenstates of both the local and the interaction 
terms. Then, for all eigenvalues £ Q (#o) associated to the eigenstates |a), we have that: 
e a (9o) — s" (9q) + ef, where e" (9q) and ef are, respectively, the interaction and the 
local contributions to the energy. On the other hand, when 9 passes from the value 
9q to 9a + 7T, the coupling contribution to the energy changes sign while the local 
one remains unchanged, so that e a (9 + n) = —ef{9) + ef. Therefore, any two energy 
eigenvalues e a {9) and e a (9 + n) corresponding to the same eigenstate \a), differ only in 
the sign of the interaction component. Then, selecting a regular step time modulation 
of the phase of the form 

fí(t\ - í 6 ' < í < T/2 ; 

periodically repeated for any t > T, we obtain that the contribution of the residual 
interaction Hamiltonian to the quantum state time evolution vanishes at any time 
t = mT with m arbitrary integer. Consequently, the effects of the undesired, residual 
interqubit XY interactions on the quantum register are completely eliminatcd, even 
in the presence of local, static imperfections in the couplings. Moreover, the result 
is vàlid for any kind of distribution of the noise on the couplings, be it Gaussian, or 
uniform, or any other probability distributution. The effects of the storing scheme are 
illustrated in Fig. Fig. and Fig. (JHJl. 

In Fig. we again plot the time evolution of the two-dimensional overlap 
J-d{t), for the same initial state considered in Fig. Q, but now under the action of 
Hamiltonian H to t(t) Eq. |JBJ. At striking variance with the unmodulated case reported 
in Fig. the step-phase modulated register realizes exact, time-periodic coherent 
revivals of the initial state. Moreover, Fig. 10} shows that the overall spatial diffusion 
of the state is always confined in a well defined and extremely narrow region of the 
ring. Because the topological control remains exact in the thermodynamic limit, the 
number of stored states con be arbitrarily large, possibly macroscopic. Finally, the 
evolution of this kind of "non-superpositional" initial state is not affected by the 
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Figure 4. Step-periodic time-modulation of the phase. Eq. 0: contour plot 
showing the evolution of the bidimensional overlap ^(í), for the same initial 
state considered in Fig. 0, as a function of the distance d from site (y-axis) 
and of the dimensionless time Aí in atòmic units h = 1 (x-axis). Here AT = 2ir 
and 9 = tt/2. The value of the overlap increases from white (^(í) = 0) to black 
Fd(t) = !)■ 




Figure 5. Step-periodic time-modulation of the phase, Eq. JJJ: contour plot 
showing the evolution of the two-dimensional overlap Tdif), for the initial state 
Ix) = (l^i) + l 1 I'-i})/v / 2, as a function of the distance d from site (y-axis) 
and of the dimensionless time Aí (x-axis). Here AT = 2tt. In graph A), 6 flips 
periodically between —tt/2 and n/2. In graph B) it flips periodically between 
and 7r. The value of the overlap increases from white ÇF^(t) = 0) to black 
ÇF d {t) = 0.5). 



choice of the initial phase 9. However, for the purposes of quantum computation, the 
most interesting and desirable scope is obviously the storage of superposition states. 
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Figure 6. Step-periodic time-modulation of the phase Eq. i71 : One-dimensional 
fidelity (t) at fixed lattice site d = as a function of the dimensionless time Aí 
(x-axis) for different constant vàlues 8o of the phase and different initially stored 
states in the case of AT = 2n. The sòlid line represents the evolution of the 
fidelity when the initial state is the single-magnon excitation l^o) and the phase 
is 6o = 7r/2. The dotted and the dashed lines both represent the evolutions of the 
initial supcrposition states |x) = (l^i) + I W - 1 } ) / %/2 . The dotted line is for the 
casc 9o = 0; the dashed line for 8o = n/2. Notice that, according to the assigned 
value of the phase, the speed of revival gain for the superpositional state can be 
larger or smaller than that of the nonsuperpositional state. 

Remarkably, if the initial state \x) is a superposition of the form 

M 

ix) = ' 

i=l 

i.e. such that the excitations are distributed among different sites, the step- 
phasc control again yields perfect periòdic state reconstruction with extremely 
limited intermediate spatial spread, although some of the details of the evolution 
at intermediate times can change significantly. As an example, in Fig. (J5J we plot the 
overlap Td(t) driven by the Hamiltonian Eq. (jüj), starting from the initial superposition 
state |x) = (|^'i) + | x I'-i))/v / 2, and for different vàlues of the phase 9. These results can 
be further generalized to any initial state \\) with an arbitrary number of excitations 
or flipped spins (arbitrarily many magnons). In this case, one finds the same coherent 
time-periodic revival of the state as in the one-magnon case, while the spatial spread 
becomes a function of the size of the region along which the initial state is extended, 
but remains in any case finite and limited. 

To better compare the different situations that can be realized using the step- 
modulated or the unmodulatcd phase, we show in Fig. © the behaviour of the 
one-dimensional fidelity at fixed lattice site d = 0, in the presence of periòdic step 
modulation of the phase for different initial states and different constant vàlues 6q of 
the phase. It is remarkable to notice that, depending on the different constant vàlues 
of the phase that one chooses as reference vàlues to operate with the control scheme, 
the approach to the full revival for the superpositional states can be larger or smaller 
than that of the non-superpositional states. An important point is to identify all the 
possible classes of states that can be perfectly stored using our method. Obviously, 
the answer to this question involves not only the interaction Hamiltonian, but the 
local, diagonal term Hq as well. If the period T of the step-phase control is gauged so 
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that 

BT = 2/tt 

with l integer, then any initial quantum state is reconstructed exactly at all times t 
integer múltiples of T. If instead BT ^ 2lir, perfect quantum state storage is still 
achieved in the subspace of all states that are linear combinations of local states with 
equal magnetization. To gain further understanding of this result, let us consider 
generalizations of the locally excited states \^d) by introducing the n-times locally 
excited states 

= n i w n i i*> ■ ( 8 ) 

j—d\,... : d n i^d\ : ... : d n 

We thus wish to analyze the case in which the state to be stored is a superposition 
of local states with different vàlues of the magnetization, and either integer or real 
vàlues of the ratio B/X. To this end, in Fig. J7J) we plot the one-dimensional hdelity 
JF (í) at fixed lattice site d — for an initially stored state of the form 

|X) = -(V2/3)|*ao) + (^ï/3)|* 72 ) + (V2/V3)|*$) , 

where the first two terms in the superposition are one-magnon states with flipped 
spin, respectively, at site d = 20 and d = 72, and the last term is the superposition 
of two one-magnon states with the excitation placed, respectively, on site d = and 
site d = 5. The initial state is thus chosen to be a superposition of states belonging 
to different subspaces of Hilbert space. In the present case, the subspace of single 
excitations and the subspace of states of the form Eq. (jHJ). As already discussed, such 
a genèric initial state achieves perfect periòdic reconstruction in the case of an integer 
valuc of the ratio B/X. In the case of noninteger ratios B/X, the state undergoes 
quasi-perfect periòdic revivals that after few cycles begin to deteriorate. Adding more 
terms to the initial superposition greatly increases the number of cycles with quasi- 
perfect reconstruction. For this kind of initial states the analysis cannot be performed 
analytically in all details, as previously done for states belonging to the same subspace. 
However, numèrics can be easily implemented for spin rings of finite, but large size. 
In the instance considered, we solve numerically the unitary evolution of J-o(t) for a 
closed chain of 90 sites. Considering larger rings yields qualitatively identical results. 



3. Other sources of noise 

Concerning the issue of practical implementations, it is crucial to verify that the 
storing scheme does not depend critically on a perfect realization of the step-phase 
modulation. To analyze this dependence, in Fig. (JHJ) we show the evolution of the 
fidclity Tç){t) — | (WolV'(í)) | 2 for the same initial one-magnon state |\&o) previously 
studied, as a function of the number of periods, when the step-periodic, time- 
modulated phase 0(t) Eq. Q is approximated by its Fourier decompositions, truncated 
at various finite orders. Remar kably, even when considering only the first 100 
harmonies, the hdelity remains close to the ideal limit J-o(t) = 1 for very long times. 
In Fig. Q we compare this behaviour with that of a superpositional initial state. 
Remarkably, the hdelity of initial superpositions, although wildly oscillatory, is always 
larger than that of non-superpositional initial states. These results prové the stability 
of the finite-harmonic approximation and, as a consequence, the robustness of the 
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Figure 7. Stcp-pcriodic time-modulation of the phase Eq. 0: One-dimensional 
fidelity .Fo(t) at fixed lattice site d = as a function of the dimensionless time 
Aí (x-axis) in a ring of 90 qubits for the initial state |x) = -(\/2/3)|* 2 o> + 
(•v/ï/3)|*T2> + (V2/V3)\^l) and for different vàlues of the ratio B/X. Plot A): 
B/X = 2, BT = 4vr. Plot B): B/X = 1.9, BT = 3.8tt. Notice the perfect periòdic 
state reconstruction in plot A) (integer value of the ratio B/X), while a small 
deviation from it causes a very slow but progressive corruption of the quantum 
state revivals, as shown in plot B). In both cases the period is fixed at the value 
AT = 2tt. 



1 . . . .. . 



0.8- 



To \ X. 

0.4- \ N. 

\ 

0.2 - \ 

v ^- 



20 40 60 80 100 

AT 



Figure 8. Fidelity J^o (í) on site i = as a function of the number of 
periods AT for the initial one-magnon state | Vfo } when the step-periodic time- 
modulated phase 8(t) is replaced with its finite-harmonic Fourier approximations, 
in increasing order. Dashed line: first 5 harmonies; dot-dashed line: first 13; 
dot-dot-dashed line: first 25; dotted line: first 50; sòlid line: first 100 harmonies. 



storing protocol against imperfections in the external control of the phase. Another 
point that needs to be addressed is the effect on the step-phase control of static, 
site-dependent imperfections in the magnètic field B that enters in the Hamiltonian 
Ho. Remarkably, it turns out that the global topological control on the off-diagonal 
Hamiltonian terms allows partial control on the diagonal terms as well. Comparing 
the time evolutions with and without step-phase modulation, in the presence of 
imperfections Si in B, with Gaussian distribution of not too broad half-width as, 
one finds that in the latter case quantum diffusion grows indefinitely (as expected), 
while in the former case it remains limited and the fidelity undergoes revivals that 
can reach unity in specific cases. Moreover, the effects of the XY residual interactions 
continue to be completely suppressed |23j . 
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Figure 9. Fidelity !Fa{t) on s ite i = as a function of the number of periods 
AT for the initial one-magnon state l^o) (black curves) versus the fidelity of the 
initial state \x) = (l^l) + |í'o))/v2 (grey curves), when the step-periodic time- 
modulated phase 9(t) is replaced with its finite-harmonic Fourier approximations, 
in increasing order. Sòlid line: first 5 harmonies; dashed line: first 25; dotted line: 
first 100 harmonies. 



4. Conclusions and outlook 

In Conclusion, we have introduced a scheme for the storage of information in a 
quantum register by a global, topological quantum control that realizes periòdic, 
perfect state reconstruction in periòdic qubit chains (rings). This kind of quantum 
register controlled by step-phase modulations is able to fully cancel the effects of 
environmental (residual) interactions of the XY type, even in the presence of local, 
static imperfections in the inter-qubit couplings. Moreover, the scheme is robust even 
in the presence of other sources of noise, such as phase modulations of finite precision 
and local static noise on the computational Hamiltonian [231 ■ The study of the effects 
of dynamic imperfections is under way |24|. However, it is already clear at this stage 
that, at least either in the ultrafast or in the adiabatic limit, the present storage 
scheme remains unaffected due to its global, topological nature. 
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